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It is commonly asserted that the electromagnetic current is conserved and therefore is not renor¬ 
malized. Within QED we show (a) that this statement is false, (b) how to obtain the renormalization 
of the current to all orders of perturbation theory, and (c) how to correctly define an electron number 
operator. The current mixes with the four-divergence of the electromagnetic field-strength tensor. 

The true electron number operator is the integral of the time component of the electron number 
density, but only when the current differs from the MS-renormalized current by a definite finite 
renormalization. This happens in such a way that Gauss’s law holds: the charge operator is the 
surface integral of the electric field at infinity. The theorem extends naturally to any gauge theory. 


I. INTRODUCTION 

A classic statement about the electromagnetic current 
in quantum field theory is that because it is conserved it 
needs no renormalization — see, e.g., 0 , P- 341], 0 p. 
430], s3j, Sec. 10.1], and Q. p. 16]. In fact this result is 
false, as we will explain in detail. Under renormalization, 
the current may mix with operators of equal or lower di¬ 
mension whose four-divergence vanishes identically, i.e., 
without use of the equations of motion. As observed in 
(Bi p. 162] there are no such operators in the absence 
of gauge fields, and so the nonrenormalization theorem 
holds in nongauge theories. The theorem also extends 
to pure QCD for the flavor currents, where the available 
operators all have nonzero color. 

However in QED (and therefore in the full Standard 
Model), the current does mix with the four-divergence 
of the electromagnetic field-strength tensor, d„F "U The 
mixing is associated with “penguin graph” contributions 
to matrix elements, Fig. mb). The conventional proof 
[I] of nonrenormalization ignores such graphs. Similar 
statements apply in more general theories 1 . 

In this article, we explain the necessary modifications, 
to all orders of perturbation theory, for the electron num¬ 
ber current in QED with a single lepton flavor. (The re¬ 
sults generalize readily to multiple lepton flavors.) There 
is a series of inter-related results: 

• The current needs ultraviolet (UV) renormaliza¬ 
tion, by the addition of a counterterm proportional 
to d v F 1 '^. The coefficient can be set in terms of the 
coupling and the Z 3 coefficient in the Lagrangian. 

• The renormalization is by addition of an operator 
that does not affect the Ward identities. Neverthe- 


1 Thus, for example, the validity of textbook statements of non¬ 

renormalization of currents for quark number, like that in 0, p. 

430], depends on whether the statements are taken in pure QCD 
or in the full Standard Model, which contains a massless photon 
field. 


less, after the application of equations of motion 
to the renormalized current operator in physical 
matrix elements, the renormalization is effectively 
multiplicative. 

• If MS-renormalization is used, the current has a 
nonzero anomalous dimension. 

• Renormalization also affects the normalization of 
the charge operator for electron number, unless the 
photon has a nonzero mass. The main part of 
this result was first found by Lurie |J p. 371, Eq. 
((8(35))]. 

• There is a unique finite correction to the countert¬ 
erm’s coefficient that both removes the anomalous 
dimension and produces a correct electron number 
operator that is the integral of an electron number 
density. In effect we have to subtract contributions 
associated with vacuum polarization. 

• Thus the Noether current for electron number is 
not the correct operator to define electron num¬ 
ber. A notable illustration is in multi-flavor QED, 
where, for example, the total muon number of a 
one-electron state is nonzero, as computed from the 
standard Noether current. 

• The counterterm restores the Gauss’s law relation 
between total charge and the flux of the electric 
field at infinity. 

• Although the effect on the total number operator, 
as opposed to the local current, only occurs when 
the photon mass is exactly zero, it depends neither 
on infrared (IR) divergences as such nor on UV di¬ 
vergences. It occurs even in a space-time dimension 
greater than 4, where there are no soft divergences 
in the scattering matrix, and with a spatial-lattice 
cutoff, which removes all UV problems. 

• A simple standard argument involving equal-time 
canonical anticommutators, appears to show that 
the electron number of states is not renormalized, 
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(b) 
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FIG. 1: (a) Graphs that are one-particle-irreducible in the current channel for insertion of a current vertex in a Green function 
or matrix element; the standard nonrenormalization argument applies only to these, (b) These graphs, one-particle-reducible in 
the current channel, also contribute to matrix elements of the current and to its renormalization. The two subgraphs that are 
cross hatched are irreducible in the photon line, while the other subgraph gives the full propagator corrections to the photon 
propagator, (c) Counterterm to (b). The filled square corresponds to an operator proportional to d„F" IL . 


contrary to reality. We resolve this paradox, which 
was first found by Weeks Q (see appendix). 

The possibility that problems exist can be motivated 
by the standard proof that the electron-number operator 
is time-independent: 



where the surface term is usually dropped. But mass¬ 
lessness of the gauge field allows a nonzero surface term. 
As we will see in detail, such a term does in fact arise, 
from penguin graphs, Fig. mb). Although the left-hand 
bubble in Fig. mb) vanishes quadratically as the external 
momentum q at the current vertex goes to zero, there is 
a pole at q I 2 = 0 in the photon propagator. In addition, 
the penguin graphs need nontrivial UV renormalization, 
with the counterterm graphs shown in Fig. |T}c) 

Although our results are quite elementary, and should 
be well-known, the only textbook account we have found 
is by Lurie p. 371]. Since current operators are aux¬ 
iliary operators used to analyze a theory, the nontrivial 
renormalization of the current does not have direct ef¬ 
fects on the scattering matrix and cross sections. Weeks’s 
paradox does symptomize some deep complications in the 
correct definition of states in a gauge theory. 

However, many practical calculations use the operator 
product expansion and factorization methods. Then ma¬ 
trix elements of currents and related operators are used, 
so that incorrect theoretical expectations for anomalous 
dimensions will affect predictions. Beneke and Neubert 
§3] indeed found closely related effects in their analysis of 


B decays, and remarked (without further comment) that 
quark number currents have nonzero anomalous dimen¬ 
sions in the presence of electromagnetic interactions. 

II. GENERAL ARGUMENT 

We use the standard gauge-fixed Lagrangian density 
£ = ■ d + eo7 • — moj ip^ 

-\{ F ^T-Y( id ' Ar (2) 

= Z 21 P (*7 • d + e/f 7 • A — too) ip 

(3) 

Our conventions are that a superscript (or subscript 
in eo and mo) denotes bare quantities, dimensional reg¬ 
ularization is in 4 — 2e dimensions, and the 7 -matrices in 
n-dimensions are normalized to obey Ti'('y ll r y u ) = 4<? M „. 
Note that the bare and renormalized couplings obey 
eo = e/A-Zg ' , with being the usual unit of mass, and 
that the gauge fixing term has no counterterm. We will 
use MS-renormalization for the interactions, although 
this will not be essential: one formula for the current will 
be in an explicitly renormalization-scheme-independent 
form. The vacuum matrix element of a time-ordered 
product of operators is denoted as j ... j ^. (The time- 

ordered product is actually the covariant T* product, as 
naturally arises with Feynman graphs or from a func¬ 
tional integral solution of the theory.) 


I 

The standard Noether current for electron number is = Z 2 'tpj tl ip, and it obeys a Ward identity: 
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Because this formula is homogeneous in the fields, it ap¬ 
plies equally whether the fields used in the Green func¬ 
tion with the Noether current are bare or renormalized 
fields; but the normalization of the Noether current itself 
is fixed. Finiteness of the right-hand side might lead one 
to conclude incorrectly that Green functions (and hence 
matrix elements) of are finite, so that the current is 
not renormalized. Certainly an extra multiplicative fac¬ 
tor is excluded [Eil, and a closely related argument 
proves that the sum of nonpenguin graphs Fig. IH a ) is 
finite: the factor Zn in takes care of divergences 

in these graphs. 

But it is possible to have a counterterm proportional 
to an operator whose four-divergence vanishes identically 
(i.e., without the use of the equations of motion). If the 
equations of motion were needed to prove the vanishing 
of the four-divergence of a counterterm, extra divergent 
terms would appear on the right-hand side of the Ward 
identity. 

In QED there is available a possible counterterm oper¬ 
ator with the appropriate dimension and symmetry prop¬ 
erties: d v F v ^. The associated divergence is in the left- 
hand bubble of Fig. mb), with its vertex for the Noether 
current. It is evidently closely related to vacuum polar¬ 
ization and hence to the wave function renormalization 
factor Zg for the photon field. The counterterm operator 
appears in graphs of the form of Fig.©c). Its coefficient 
can be determined by the operator equation of motion 
for the renormalized photon field: 

0= JTTT = e ^N + ^ , ‘ + \d*d-A, (5) 

where S = f d 4 ~ 2e x C is the action for the theory. Hence, 
uniqueness of expansions in poles at e = 0 shows that the 
MS-renormalized current is 


4s = ^ + (6) 

in terms of which the photon equation of motion has a 
simple form in terms of finite operators: 

0 = e^j» m + d v F^ + h»d-A. (7) 


Now the use of the operator equations of motion in¬ 
duces extra delta function terms when applied in Green 
functions. However this is not the case in matrix ele¬ 
ments. Moreover, the Gupta-Bleuler condition gives zero 
physical matrix elements for d ■ A. So we can eliminate 
either one of the operators in Eq. © in favor of the other, 
plus a term that vanishes in physical matrix elements: 


7— 

■'MS 


d v F v » 

e/j, e 


—• A. 



1 — 7~ 1 

- Z-d^d-A. 


( 8 ) 

( 9 ) 



FIG. 2: One-loop wave-function and vertex corrections. 



FIG. 3: One-loop penguin graph contributing to the renor¬ 
malization of the vector current. The cross denotes an inser¬ 
tion of the vector current. 


Thus in physical matrix elements, the renormalized cur¬ 
rent obeys: 

r? 1 

•?MS = - U ~~T = ( in Physical ME). (10) 

C/i z^3 

The formulae with only gauge fields exhibit the finite¬ 
ness of the renormalized current, while those with the 
Noether current exhibit the nontrivial renormalization of 
the current 

One further useful result is the expression of the cur¬ 
rent in terms of bare fields: 

= ^ (0) yv (0) + "F (11) 

Mb eg 


The multiple formulae for the current lend themselves 
to different interpretations. Within Feynman graph cal¬ 
culations, the distinction between a counterterm propor¬ 
tional to d u F v,Jj and one proportional to is absolutely 
clear and unambiguous: Thus the left-hand bubble of 
Fig. mb) is renormalized using the d v F v>l operator, not 
But after applying the equations of motion, the dis¬ 
tinction is not so clear cut. Indeed, using the formula for 
the renormalized current in terms of the Noether current 
will lead us to Weeks’ paradox in Sec. El 


III. ONE-LOOP VERIFICATION 


In a one-electron state, the one-loop matrix element of 
the electron number current has contributions from wave- 
function and vertex graphs, Fig. ©I and from a penguin 
graph Fig. ©I As is well-known, the UV-infinite parts of 
the wave-function and vertex graphs cancel. Moreover, 
at zero momentum transfer in an on-shell state, the com¬ 
plete vertex and wave function graphs cancel, so that 
they make no contribution, for example, to the expecta¬ 
tion value of the electron number operator. 

However, the penguin graph also contributes. Its 1PI 
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part is 


be related to the RG of the interaction by use of Eq. CD: 


eS e )J, 6 
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where S f = (4ne TE ) e . The graph has a divergence which 
is canceled by the MS-counterterm in Eq. © with 


Z 3 = l- 


e 2 S e 

127r 2 e 


+ 0(e 4 ). 


(13) 


Here we use the usual QCD definition of MS- 
renormalization that the counterterms contain a factor 
of S € for each loop. 


IV. APPARENT ELECTRON NUMBER 
ANOMALY 


From our calculations, we can see that the penguin 
graph contribution to the matrix element of the current 
apparently changes the value of electron number, no mat¬ 
ter whether we use the Noether current itself or the MS- 
renormalized current. For the total electron number we 
need the limit as g 2 —► 0. Although the loop vanishes 
at q 2 = 0 , this is canceled by the pole in the photon 
propagator. Thus the matrix element of jiig in a single¬ 
electron state at zero momentum transfer (and e = 0 ) 
is 


= 




2 (14) 


so that the electron number of an electron is apparently 


1 + 




0 (e 4 ). 


(15) 


This is not even renormalization-group invariant. There¬ 
fore we cannot interpret the current as corresponding to 
electron number in the standard way, despite the fact 
that the current has the correct Ward identity that cor¬ 
responds to the standard commutation relations between 
the current and the fields. 

The structure of the graphs in Fig. U shows that the 
factor m is universal between different matrix elements 
of the current and charge. 

An even more dramatic and paradoxical consequence 
occurs if we add a second flavor of fermion, which we can 
call a muon. The matrix element of the muon-number 
current has a contribution from a graph of the form of 
pig. a with the electron loop replaced by a muon loop. 
Thus the muon number of the electron appears to be 
nonzero. 

It is also useful to examine the RG properties of the 
MS-renormalized current. Its anomalous dimension can 


u A 1 n u u r r, u v r n 

'V® = 2lA ~Fzr = 2ia —’ (16) 

where 7,4 is the anomalous dimension of the photon field: 


1 dZ 3 e 2 S e 


7A(e) = = 


(D(e 4 ), 


(17) 


2Z 3 d/r 12 - 7 T 2 
which is related to the QED /3-function by 
de 

=-ee +P{e) =-ee +ejAie). (18) 


The nonzero anomalous dimension of the current, in 
Eq. Ill (il l, is consistent with the renormalized equations 
of motion, Eq. ©. We see this by use of the relations 
between /3 and 7 , 4 , together with the anomalous dimen¬ 
sion of the gauge-fixing parameter, d£/d/r = — 27 , 4 : 


0 = ^( d " F ‘ / “ + e i‘‘im + \ d * d ' A ) 

= 7,, - 1//'.)';.. + A 1 ;) - .1 j (19) 


The right-hand side is proportional to the equation of 
motion, so that it vanishes. Finally, applying Eq. CE 3 
gives a renormalization group equation for the current in 
terms of itself when physical matrix elements are taken: 

= ~ 2 7a ( in physical ME). (20) 


V. REDEFINITION OF CURRENT 


Since the apparent anomaly arises from penguin-graph 
contributions only when the photon is exactly massless, 
a more satisfactory definition of the current evidently 
requires us to remove the penguin graph contributions. 
To preserve the locality of the current, we can perform an 
exact removal only at one value of q , naturally q = 0. We 
will demonstrate that the correct definition (only valid if 
the electron mass m is nonzero) is: 


= j— _- 

J •'MS 


n(0) d v F v » 

e/A 


— Jn + [^3 ^ 
= 


1 - n(o)] 


1 - 


d v F v » 

e/A 

l + n(0)' 


d v F (°) ^ 


eo 


( 21 ) 


Here n(g 2 ) is the vacuum polarization, defined as usual 
so that the renormalized photon propagator is 


i (-9nv + q^v/q 2 ) iq^quC 

q 2 [1 + n(< 7 2 )] (<? 2 ) 2 
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The first form of j M in Eq. m shows that it is a finite 
operator that obeys the standard Ward identity, Eq. iS- 
In the last form, the factor [1 +11(0)] /^3 is the inverse 
of the photon-pole residue in the propagator of the bare 
photon field. Thus the last form is written solely in terms 
of bare quantities, so that the current is RG invariant. 
Moreover we have a formula where a nontrivial correction 
is manifestly needed even if there are no UV divergences, 
as when n < 4. 

The extra photon term in the definition of the current 
depends on the dynamics of QED, so we call it the dy¬ 
namical term in the current. Its normalization is only 
known after the theory is solved, and depends on specif¬ 
ically quantum mechanical effects. 

In physical matrix elements, the equation of motion for 
the photon field gives 


_ [1 + 11 ( 0 )] d u F^ 

e/+ 

= I 1 + n (0)] 
i + n(o) 

^3 Jn 


d F vM 

^^phys 

^phys 


(all in physical ME). 


(23) 


In the last line, we use what we term a physically nor¬ 
malized field, 


iX. = [l + n(0)] 1 / 2 iW‘, (24) 


whose photon pole has unit residue. The corresponding 
coupling is 


e P hys = [l+n(0)]- 1 / 2 e M £ , (25) 


which has the conventional value, i.e., 1.6022 • • • x 
1(T 19 C, rather than the MS-renormalized coupling e 
that is appropriate in certain high-energy calculations. 
Since, with MS-renormalization, 

n(0) = -li ln ^ + O(e4) ’ (26) 


current conservation shows that only the first term in Eq. 
123 survives, so that the matrix element differs from the 
irreducible contribution by a factor [l + n(0)]/[l + II(g 2 )], 
which is unity when q 2 = 0. The standard arguments 
about nonrenormalization of the current actually apply 
only to the 1-photon-irreducible part, Fig. which is 
finite and unaffected by the addition of the dynamical 
photon term to the current. 

To see that the extra dynamical term in the current 
is due to bad behavior of the current operator at spatial 
infinity, we can use a nonzero photon mass. Then the 
factor between the total and irreducible parts of a matrix 
element of the current becomes 

q 2 -m 2 + g 2 n(0) g 2 [n(0) - n(g 2 )] 

g 2 — m 2 + g 2 II(g 2 ) g 2 -m 2 +g 2 n(g 2 )’ 1 ’ 

where the last term, proportional to 11(0) — n(g 2 ), is the 
contribution of the penguin graphs and of the dynamical 
term. Because of the photon mass in the denominator, 
both of these contributions vanish at g = 0. 

Finally, we observe that the formula JSU for the phys¬ 
ical current in terms of the electromagnetic field strength 
shows that the time-component of the current is 


e p hys 


This is the divergence of the electric field operator, with 
its standard normalization, in units of the (negative) 
charge of the electron. Integrating over all space shows 
that the electron number with this definition equals the 
value given by Gauss’s law with the standard normaliza¬ 
tion. This naturally matches with the classical limit of 
electromagnetism, for macroscopic phenomena. 

We have defined a physical charge and current, but 
the necessity of nontrivial renormalization depends on 
whether the photon is massless or not and on whether 
the current or charge is considered: 


• In the 3+1 dimensional theory, the Noether current 
always needs UV renormalization, independently of 
m 7 . 

• The counterterm never affects the Ward identities. 


we verify that the redefinition Eq. El does indeed re¬ 
move the one-loop anomaly in the electron number. 

To see that the current j^ defined in Eq. El is the 
uniquely correct current, we express an arbitrary Green 
function or matrix element of in terms of Gq j, which 
is the 1-photon-irreducible part in the current channel. 
The reducible graphs are those shown in Fig. mb) and 
(c), so that the total Green function is 


'-'tot 


1 + 11(0) , ^ q^ n(g 2 ) - n(o) 
ui 1 + n(g2 ) +Ct ui q 2 i + n(g 2 ) ‘ 1 ] 


Multiplying by g M gives just the irreducible contribution 
which verifies that the Ward identities are un¬ 
affected by the extra terms used to define the physical 
current. But when we take a physical matrix element, 


• When the photon mass is nonzero, the counterterm 
integrates to zero, so that the charge does not then 
actually need renormalization, and the integral of 
the Noether charge density gives the correct charge. 

• But when the photon mass is zero: 

— Integrating the counterterm over all space 
gives a nonzero result, so that the charge needs 
UV renormalization. 

— A particular finite part in the counterterm is 
needed to obtain the correct charge. 

— Even when the theory is regulated in the UV, 
the finite renormalization of the current and 
charge are still needed, if the the charge is to 
be correct. 
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• However, even with a nonzero photon mass, the def¬ 
inition of the physical current m is valid and gives 
the correct charge. (The counterterm integrates to 
zero.) This definition has the practical advantages 
that there are no problems in taking the limits of 
zero photon mass, integration over all space, and 
removing a UV regulator, and that changing the 
order of limits is safe. 


Evidently the need for a nontrivial renormalization of 
the total electron number (as opposed to the current) 
only arises when the photon is massless. However, it is 
to be emphasized that this is not due to actual infra-red 
divergences as normally understood. This can be seen by 
considering the theory in 5 space-time dimensions (i.e., 
4 space dimensions), with a spatial lattice. The higher 
space-time dimension is sufficient to remove all the usual 
soft divergences in the scattering matrix. The use of an 
integer dimension removes all possible artifacts associ¬ 
ated with dimensional regularization, and the use of a 
spatial lattice gives us a conventional quantum mechani¬ 
cal theory (in real time) without UV divergences. All the 
considerations leading to a physical current that is not 
equal to the Noether current still apply. 

It is sometimes claimed that the actual definition of 
QED requires us first to regulate the theory in the IR, 
say with a nonzero photon mass, and then to take the 
limit of zero photon mass. If this were so, we could avoid 
the need to redefine the charge, since in the regulated 
theory the Noether charge would give the correct answer. 
However, as a field theory, QED exists if the photon mass 
is kept at zero at all stages; this is evidenced by the fact 
that the Green functions of the theory exist directly at 
m 7 = 0. The use of an IR regulator (and inclusive cross 
sections that allow undetected soft photon emission) is 
only necessary if one uses the conventional LSZ formal¬ 
ism to compute cross sections. In any case, when the 
space-time dimension is above 4 there are no soft photon 
divergences in cross sections, but nevertheless a redefini¬ 
tion of the charge is still needed, as we have just seen. 

Our arguments generalize to more complicated theo¬ 
ries. For example, the approach applies to QED with 
more than one flavor of lepton, with suitable changes in 
the numbers of flavors in the vacuum polarization graphs. 
The approach applies not only to the electromagnetic 
current itself, but also to the individual conserved flavor- 
number currents. 

We propose an interpretation of the nontrivial renor¬ 
malization of the current as a universal effect of attempts 
to measure the current. The quantum mechanical current 
creates polarization of the vacuum. Although this effect 
is proportional to 1/r 2 at a distance r from a source, it 
has to be integrated over a sphere of surface area propor¬ 
tional to r , so that the effect is nonvanishing as r —> oo. 
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APPENDIX A: WEEKS’ PARADOX 

1. Formulation 

Conventional expectations for the interpretation of a 
charge operator, and, in particular, for its value in a state 
of definite electron number arise quite generally from the 
commutation relations between the Noether current and 
bare elementary fields. These commutation relations are 
entirely unaffected by UV renormalization and are fun¬ 
damental to the definition of a quantum field theory. In¬ 
spired by the paper of Weeks [jJ, we now present a para¬ 
dox. 

The paradox and its resolution turn on the issue of the 
correct way to implement concretely the concept of the 
electron number of a state. There are several possible 
ways of implementing the concept. In scattering theory, 
we can use in- or out-states as basis states; the labeling 
of these states in the usual way gives unambiguous values 
for electron number. But we could also define electron 
number in terms of eigenvalues of a suitable operator for 
conserved charge. Yet another possible definition would 
use the expectation value of the charge in a state. In 
the absence gauge fields, the different definitions can be 
proved to agree. But, as we will see, in QED an attempt 
to make such a proof runs into problems. The problems 
are related to, but distinct from, the issues treated in the 
body of this paper. 

The paradox is obtained as follows: 

1. Let |0) be the (true) vacuum state and let |1) be any 
state of electron number unity. While we may use 
an on-shell one-electron state of definite momentum 
for 11), it is safer to use a normalizable wave-packet 
state. Here we chose the states to be normal in- 
or out-states, with electron number being defined 
with respect to the obvious labeling of the states in 
terms of particle content. 

2. We assume, as is natural, that (1|^|0) ^ 0, and that 
|0) and |1) are eigenstates of the Noether charge 
Q n- Of course, we have shown in the body of 
this paper that there is an unexpectedly nontriv¬ 
ial renormalization of the charge, which leads to 
the expectation that the eigenvalue disagrees with 
naive expectations. But this finding by itself does 
not impinge on the presumption that we have an 
eigenstate. 
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3. We (and Lurie @) have calculated that the expec¬ 
tation value value of Q n in the state |1) is Z$ . Here 
Z 3 is the renormalization factor for the photon field 
in an on-shell field. It equals Z^/[\ + n( 0 )] in our 
earlier notation. An explicit one-loop calculation 
shows that Z 3 ^ 1. Of course, on an eigenstate the 
expectation value of an operator equals the eigen¬ 
value. 



FIG. 4: Matrix element on the right of Eq. used in 

lowest-order BJL computation of the commutator. The cir¬ 
cled cross is the vertex for the operator and the simple 
cross is the vertex for tp. The external line of momentum p is 
for the state |1). 


4. From the canonical equal-time anticommutation re¬ 
lations for the fermion field we find 

(l|^(t,x)|0) = (1| [Qn, x)] |0> 

= <1|Qn#)>-<1|^Qn|0) 

= (z 3 -0) (l|#,x)|0) (Al) 

Since Z 3 ^ 1, there is a contradiction, i.e., there is at 
least one mistake in the calculation or the assumptions. 

Weeks’ 0 derivation of the paradox was in Coulomb 
gauge, but the main ideas are the same. The basic 
conflict is between the canonical commutation relation 
[<2n> $] = 1 • and the value Z 3 1 for Qn in the state 
| 1 >- ' 


2. Resolution 


The resolution of the paradox is that physical states 
are not exactly eigenstates of the Noether charge oper¬ 
ator Q n , contrary to the (unproved) assumption made 
at step 2 of the derivation of the paradox. For example, 
as we will see, when the charge density operator j’fc is 
applied to the vacuum and then integrated over position, 
there is a certain contribution that vanishes in the limit 
of infinite volume (or equivalently when the momentum 
transfer q goes to zero). However, matrix elements of a 
field operator with this part of the state have a divergence 
in the same limit. Moreover this anomalous contribution 
is associated with unphysical parts of the theory: the un¬ 
physical states in Feynman gauge, or the instantaneous 
Coulomb potential in Coulomb gauge. Thus it is not 
easily visible when one restricts attention to physical ma¬ 
trix elements; in particular, the anomalous contribution 
does not affect the calculation of the expectation value 
of the charge in physical state. In addition, in Coulomb 
gauge the equal-time commutators and anticommutators 
are modified by interactions. 

We will illustrate these issues by calculations of the 
commutators in two ways. One is the Bjorken-Johnson- 
Low (BJL) method 9], which is applied to time-ordered 
products of the current and fields, while the second 
method is a direct calculation with the intermediate 
states used in the commutator — see Eq. 63 below. 



FIG. 5: Graph for the one-loop vacuum-polarization contri¬ 
bution to the BJL computation. 


a. Resolution in Feynman gauge 


Consider first the commutator in covariant gauge be¬ 
tween the Noether current and the electron field. 

( 1 l[jN(*>y)>'*/’( i - x )]|o) = 

-^<lHr}<Yb». (A2) 

Y 

The sums over intermediate states are over all states, not 
just over representatives of the physical subspace. 

The BJL method gives this matrix element of the com¬ 
mutator from the Fourier transform of the corresponding 
matrix element of the time-ordered product, by taking 
the limit as q° —► 00 of iq° times the matrix element. 
Here q M is the external momentum flowing into the ver¬ 
tex for Thus 

(ilbN^w),^,*)]^) = 

lim iq° f dt e iq0t (l\T j^(t, y), ip(t, x)|0). (A3) 

9°->oo J 

It is readily verified that the lowest-order calculation 
from the graph of Fig.Q]gives the expected commutator: 


lim 

q °—»oo 


iq°u 7 0 


i(j> — fa + m) 
(p — q ) 2 — m 2 


= u, 


(A4) 


where u is the Dirac wave function for the on-shell one- 
electron state 11), of momentum p. 

The one-loop penguin graph contribution to Eq. lA3l) 
is from the graph of Fig. [3] Any correction factor to the 
LO result is from the large q° behavior of the vacuum 
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FIG. 6: Graphs for the one-loop vacuum-polarization con¬ 
tribution to the BJL computation, with cuts on the photon 
line. 



FIG. 7: Graphs for the one-loop vacuum-polarization contri¬ 
bution to the BJL computation, with cuts on an electron- 
positron pair. 


polarization graph times the photon propagator: 

+ 9 V)n(<z 2 ) = (V° + n ( q 2 ). 

(A5) 


The [i = 0 case is 


(A6) 

which is suppressed at large q° by two powers of q°. The 
spatial part is 


^V), (A7) 

q 2 

again suppressed, but by one power of q°. 

So the canonical commutator is unchanged. But notice 
how the q^q° part of the vacuum polarization is essential 
to the suppression. This contrasts with the fact that this 
term does not contribute to physical matrix elements of 
the current. We may therefore anticipate that contribu¬ 
tions from unphysical intermediate states in Eq. are 
essential to get the correct result. 

This can be seen explicitly from the graphs of Fig. 0 
where the cuts correspond to the cases that |X) is an 
out-state of one electron and one photon, and that |Y) is 
a state of one photon. The vacuum polarization factor is 


K-g^q 2 + <zV) n(g 2 ) = <|q|n(o), (as) 

given that the photon is on-shell, q 2 = 0. This means 
that in the first graph, the sole contribution is from a pho¬ 
ton with a scalar polarization, i.e., a polarization vector 
proportional to its lightlike momentum. This is a state 
of zero norm: the current applied to the vacuum has not 
created a genuine photon. Moreover, in the relevant part 
of IA2I) . the vacuum polarization factor (IA8I) vanishes 


quadratically when q —> 0, which is completely compat¬ 
ible with the natural expectation that the vacuum has 
zero charge, so that it is annihilated by the operator Q n- 

In the commutator, this factor contributes to the 
matrix element (y|j^|0). However it is multiplied by 
<l|^|y), and then we get a nonzero result from the fol¬ 
lowing dependence on q: 

• The already-derived factor of q 2 . 

• 1/1q| in the photon relativistic phase space. 

• The IR singularity oc l/|q| in the fermion propaga¬ 
tor in (l\ip\Y). 

The nonzero result at q = 0 is independent of the space- 
time dimension. Thus it does not depend on the existence 
of IR divergences that violate the standard principles for 
asymptotic scattering states, but only on the massless¬ 
ness of the photon. 

To summarize: Although the charge operator applied 
to the vacuum, for example, does give zero, this is com¬ 
pensated by an infinity in the matrix element with the 
field operator. Moreover the offending intermediate state 
is of zero norm, before taking the q —> 0 limit, so the state 
is equivalent to the zero state. Similar remarks apply to 
the one-electron-one-photon intermediate state from the 
other diagram in Fig. El 

A related issue is about the value of the commutator 
of the physical charge. In physical matrix elements, the 
two charges differ by a factor Z%, and it is tempting to 
say that the commutators also differ by this factor. But 
consider the expression of the physical current in terms of 
the Noether current and d ■ A. This is obtained by appli¬ 
cation of the equations of motion Eq. 0 to the physical 
current Eq. <EUk but without the physical state condition 
that we used in Eq. lf23l) : 

r= j- 3 fK+ %i^ dfdA <A9) 














9 


The physical electron number operator is defined as an 
integral over j°. so that 

Qphys = i Qn + Z \ [ d 3 4 x 8 ° {d ■ A). (A10) 

Zj, J 

In physical matrix elements we can ignore the second 
term, to get the usual ratio of Zj, between the physi¬ 
cal and Noether charges. But in commutators with ele¬ 
mentary fields, we must include the second term. Since 
this second term involves a double time-derivative of the 
gauge field, the equations of motion must be applied to 
write it in terms of elementary fields and their canonical 
momenta, if we are to use canonical commutation and 
anticommutation relations. This gives extra terms, so 
the commutators of Qphys and Qn with elementary fields 
are not simply related by Z$. 

The appropriate formula to use is, in fact, either of the 
two last lines of Eq. m • For the /i = 0 component we 
have 


^?phys — Q N + 


1 - Zf 1 


= Qn ~ 1 


-l 


d 3 x 


o F (°) uO 


eo 


eo 

V 2 A° + 4-V • A (0) 
at 


(All) 


terms in the propagator that are exactly zero because the 
vacuum polarization is transverse. It is easily checked 
that this formula equals — g^A(q 2 ). As q° —* oo, this is 
missing the power suppression that we have in covariant 
gauge. If the theory is regulated in the UV, e.g., by using 
e > 0, then the vacuum polarization does go to zero, 
so that the loop correction to the standard commutator 
vanishes. But after the regulator is removed, the vacuum 
polarization grows logarithmically. 

Evidently there is an anomaly in the commutator of 
the Noether charge density with the field ip. Because 
constrained quantization is used, it is not automatically 
true that standard canonical commutation and anticom¬ 
mutation relations can be maintained. But this result 
depends on the space-time dimension. 

Interesting dimension-independent pathologies also 
arise in the calculation of the commutator with inter¬ 
mediate states, but from electron-positron terms instead 
of photon terms. In the photon term, we have the same 
vacuum polarization factor Eq. (IA8I) as before, but now 
q M times the numerator on the on-shell photon line gives 
exactly zero. 

Instead there are extra contributions from intermedi¬ 
ate states with electron-positron pairs, Fig. 0 The cut 
vacuum polarization graph has the transverse form 


i.e., the Noether charge density, with unit normalization, 
plus a term that involves at most a first-order time deriva¬ 
tive, so that it commutes at equal time with the fermion 
fields, to give [Qphysj ■*/>] = [Qnj 1 /’] = 4>- Thus, despite 
the relative factor Z% between the physical and Noether 
charges as projected onto physical states, the commuta¬ 
tors with the fermion do not have this factor. 


b. Resolution in Coulomb gauge 


In the Coulomb gauge version of the BJL calculation 
from the graph of Fig. |3 the vacuum polarization graph 
times the photon propagator is now 


+ q^v^q 0 / q 2 ) 


*(-^V+ 9V)n( 9 2 ), 


(A12) 


where v M = ( 1 , 0 ) is the unit vector defining the rest 
frame for the Coulomb gauge, and we have dropped those 


(-^Y + q v q 0 )U c (q 2 ) = (q 2 , q 9 0 ) , (A13) 


where q 2 is no longer zero, but is bounded below by 4m 2 . 
The spatial part of this vector gives zero when multi¬ 
plied into the photon propagator. However, the time 
component multiplies a factor j/q 2 for the instantaneous 
Coulomb potential, so that we get a nonzero limit at 
q = 0. Of course, the energy q° is nonzero, but this 
is allowed, since the charge operator is computed from 
the charge density at a fixed time, so that in momentum 
space there is an integral over all q°. 

The result of a nonzero limit at q = 0 is dimension 
independent, and again indicates that states of definite 
particle number are not actual eigenstates of the Noether 
charge. But the source of the anomaly has changed. In 
Feynman gauge the anomalous contribution was associ¬ 
ated with unphysical states. But in Coulomb gauge, it is 
associated with the unphysically instantaneous potential 
in the photon propagator. 
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